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1 Introduction 

Ball mills have been deployed as a simple top-down approach for comminution processes up 
to nanometer scale. Although its practicality, ball mills involve many parameters related 
to both external and internal dynamics of the equipments. All of them unfortunately lead 
to severe uncertainties and become some major obstacles to perform an efficient milling. 

Meanwhile, more theoretical approaches based on mathematical description are ex- 
pected to overcome such problems. Many efforts have been done to develop modelization 
of ball mill equipments [1, 2, 3, 4, 5]. Through the modeling approach and its subsequent 
simulation, one expects to be able to obtain prior information and constraints to optimize 
experimental strategy. 

In general, the experimental works using ball mill equipments face the following prob- 
lems : 

• Choosing an appropriate parameter set related to the milled materials for particular 
case of experiments. It could be the ball size, the initial size of milled powders, the 
number of ball, the initial quantity of powders, material characteristics of ball and 
also powders, and so forth. 

• The milling time for particular characteristics and quantities of balls and powders. 

• The design and geometrical motions of vial itself. Although there are various types 
of ball mills, most of them have not been developed based on prior comprehensive 
simulations. 

• On the other hand, the experimental measurements on the internal dynamics of vial 
are almost impossible. 

The question is then how to overcome these problems ? Is there any smart solution for 
these obstacles ? 

In our recent works, a novel model to describe the internal dynamics and to relate it 
with surrounding environment has been proposed [6, 7]. In contrast to the semi-empirical 
approaches, the model does not require prior experimental data or simulation results to fit 
the parameters [8, 9, 10, 11], nor huge computational power as in some more empirical ap- 
proaches [12, 13, 14]. Actually our model combines the deterministic approach for milling 
bodies motion, and the statistical approach to relate them with external macroscopic 
physical observables, in particular system temperature. 

In this paper, however the focus is put on investigating the behaviors of normalized 
macroscopic pressure, P/Pq. The observable is important and accessible in most of real 
experiments. Particular interest is investigating its dependencies on the system tempera- 
ture and the evolution milled powder mass as well. 

The paper is organized as follows. First, after this introduction the model is briefiy 
explained. Before summarizing the results, numerical analysis and simulation for the 
normalized pressure in term of system temperature are discussed. 

2 Model and simulation 

Rather solving a set of equation of motions (EOMs) governing the whole dynamics as 
always done in conventional approaches, the dynamics is described in a hamiltonian in- 



volving all considerable potentials in the system and its surrounding environment. 

In the model, the dynamics of each 'matter' in the system, i.e. balls and powders 
inside the vial, are described by a hamiltonian i/m(^, t). The index m denotes the powder 
(p) or ball (b) and r = {x,y,z). The hamiltonian contains some terms representing all 
relevant interactions working on the matters inside the system as follow [7], 

Hm = Hq + v/m-m + Kn-v + Kn-m' + Kxt 5 (1) 

with V denotes the vial, while Hq is the free matter hamiltonian containing the kinetic 
term, 

i=l 

where n^ is the matter number, m^ and iv are the matter mass and momentum respec- 
tively. Throughout the paper we assume that the mass or size evolution of matters is 
uniform for the same matters. 

The matter self- interaction Vm_m, the matter-vial interaction Vm_v and the inter- 
actions between different matters may be induced by, for instance, the impact (\/^™p) 
potential, 

K^'l^ t) = - E E / """^ ^^ d (W)., n . (^3)^ . , (3) 

with m, m' : v, p, b and n is the unit normal vector. The potential should in fact represent 
the whole classical dynamics among the matters, i.e. the impact forces among balls and 
powders. Here the impact force is dominated by its normal component [14], 
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Here, T^„^/ is the Young modulus, f mm' represents the Poisson ratio of the sphere material, 

-^mm' — {RmRra') / {Rra + Rra') IS the cffcctive radiuS, whilc .^mm' = Rm + Rm' — l^m" "^m'l IS 

the displacement with R^ is the radius of interacting matter. A is a dissipative parameter 
given in [15, 16, 17] containing the viscous constant r]^. 

In fact, there are another potentials like Coulomb and gravitational potentials which 
may influence on the system. However, those contributions should be considerably tiny 
due to its neutral charges and tiny masses as well [7]. 

Incorporating the effect of external electromagnetic field surrounding the system for 
charged matters shifts the kinetic term in Eq. 2 to be [7], 

Ho — >Ho+EM = - '^\iPni)i- QmA +n^Q^(t), (5) 

1=1 

with electromagnetic (scalar and vector) potential (0, A) and the matter charge Qm- 

From now, let us focus only on the dynamics of powders which is our main interest in 
the sense of comminution process. From Eqs. (1), (2) and (3), the total hamiltonian for 
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for Qp 7^ 0. The last two potentials represent the total impact potential among powders; 
powders and vial; powders and balls respectively. Obviously we do not need to take into 
account the ball self-interaction Vj^™^ '^*-'^ ball-vial interaction V^^. This is actually the 
advantage of using hamiltonian method. 

Relating a hamiltonian with macroscopic physical observables can be realized through 
the partition function known in statistical mechanics, 
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for a canonical ensemble of matter m governed by a particular hamiltonian H^. Here, 
(5 = l/(kBT) with kB and T are the Boltzman constant and absolute temperature. Further 
one can calculate, for instance the normalized pressure as, 

PL='^. (8) 



Integrating out the time component at finite temperature, one obtains [7], 
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Eq. (9) provides a general behavior for temperature-dependent pressure in the model, 
while the geometrical structure and motion of vial is absorbed in the function J-'. 

This is the final result which is ready to be evaluated further using numerical ap- 
proaches like Monte Carlo. 
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Figure 1: The normalized macroscopic pressure, P/Pq, 
typical order values of J^ = 0{lO-^^) ~ 0{lO~^^). 



as a function of temperature for 



3 Results and summary 

In this paper, the integral in Eq. (10) is performed using Monte Carlo technique. Sim- 
ulation on the system temperature dependencies of normalized pressure is shown in Fig. 
1. Secondly its dependencies on matter mass is given in Fig. 2. 

The simulation is done for steel material in a typical geometry of spex mixer / mill 
with resolutions on 3-dimensional space ~ 10^ x 10^ x 10^, vial length / = 50 mm, vial 
radius r^ = 10 mm, shaft-arm length L = 200 mm and ball radius R^ = 5 mm. 

From the figures one can conclude that high efficiency is expected at low system tem- 
perature region as can be understood from natural sense since the higher is temperature, 
the agglomeration phenomena occurs which prevents further comminution processes. In 
another words, the higher is temperature, the pressure is getting down which lead to lower 
impact energy. Therefore these results suggest that in general one should keep optimized 
working temperature as long as possible to keep high energy milling. 
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Figure 2: The normalized macroscopic pressure, P/Pq, as a function of powder mass for 
J-" ~ O(10~^^) at typical room temperature. 
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